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Abstract 

We analyze the twist-4 contributions to Bjorken and Ellis- Jaffe sum rules 
for spin-dependent structure function gi(x, Q 2 ). We investigate the anoma- 
lous dimensions of the twist-4 operators which determine the logarithmic 
correction to the 1/Q 2 behavior of the twist-4 contribution by evaluating 
off-shell Green's functions in both flavor non-singlet and singlet case. It 
is shown that the EOM operators play an important role to extract the 
anomalous dimensions of physical operators. The calculations to solve the 
operator mixing of higher-twist operators are given in detail. 
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1 Introduction 



The polarized structure functions g\ and (72 for nucleons are measured by recent exper- 
iments at CERN [I], ^| and SLAC f||, |]]. These functions provide us the non-trivial spin 
structures of nucleons. Especially g\ function has direct partonic interpretations, and it 
turned out that very relativistic pictures hold for nucleon spin |J. In the framework of 
the operator product expansion and the renormalization group method based on QCD, 
we can derive sum rules for n-th moments of gi' n (x,Q 2 ), where the first moments of 
gi' n (x, Q 2 ) are given by the Ellis- Jaffe sum rules and that of g v ^ n {x, Q 2 ) leads to 
the Bjorken sum rule 0. In the deep inelastic scattering, the perturbative QCD has 
been tested so far for the effects of the leading twist operators, namely twist-2 opera- 
tors, for which the QCD parton picture holds. Now, the twist-4 operators give rise to 
0(1/Q 2 ) corrections to the first moment of gi(x,Q 2 ), which may have some contribu- 
tion in the region of Q 2 of the recent experiments. Their contributions correspond to 
the correlations between quarks and gluons. 

In this paper we investigate the renormalization of the twist-4 operators, which are 
relevant for the first moment of the nucleon spin structure functions g\(x,Q 2 ), and 
obtain their anomalous dimensions. From this calculations we determine the logarith- 
mic correction to the 1/Q 2 behavior of the twist-4 operator's contribution in the first 
moment of g\(x,Q 2 ). We study the renormalization of the composite operators by 
evaluating Green's functions taking the external lines off-shell so that we can avoid 
the subtle problems of IR divergence. The general feature in renormalization of higher 
twist operators by off-shell Green's functions is that there occurs the mixing among 
the physical operators and the EOM operators which are proportional to the equation 
of motion and what we call the 'BRS-exact' operators which contain ghost operators 
H, |J. Although the physical matrix elements of EOM and BRS-exact operators van- 
ish, we need to consider properly their contributions to the counterterms to extract the 
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anomalous dimensions of the physical operators. As we show later, the tree vertices of 
these unphysical operators may have components of the same tensor structures as that 
of the physical one. So we need to identify correctly the operator basis and separate the 
divergence of the radiative corrections into two parts that correspond to the physical 
operators and the unphysical operators, respectively. As for the operator mixing, some 
formal arguments can be given in the gauge theory || ||, and in this paper we confirm 
these theorems explicitly. 

In the following sections, we obtain all mixing matrix elements for flavor non-singlet 
part, and only the physical one for flavor singlet part. In both cases, it turns out that 
there exists only one physical operator 



2 Twist-4 contribution to the first moment of gi(x, Q A 

The polarized deep inelastic scattering is described by the antisymmetric part of the 
hadronic structure tensor W^ v given in terms of two spin structure functions g±(x, Q 2 ) 
and gi{x, Q 2 ) as 

= e, uXa q x j^-i-^x, Q 2 ) +(p. q8 '- q . sp °)-L-g 2 ( Xj Q2) j , (l) 

where q is the virtual photon momentum and p is the nucleon momentum, x is the 
Bjorken variable x = Q 2 /2p ■ q = Q 2 jlMv and q 2 = —Q 2 . M is the nucleon mass and 

= u(p, s)7 At 75w(p, s) is the covariant spin vector of the nucleon. 

Now the first moment of the gi(x,Q 2 ) structure functions for proton and neutron 
turns out to be up to the power correction of order l/Q 2 : 



rr(Q 2 )= f 9 r(x,Q 2 )dx 

Jo 



where gi^ n \x,Q 2 ) is the spin structure function of the proton (neutron) and the plus 
(minus) sign is for proton (neutron). On the right-hand side, a 3 = gA = Ga/Gv is the 



ratio of the axial- vector to vector coupling constants. Here we assume that the number 
of active flavors in the current Q 2 region is Nf = 3. Denoting (p, s\ipj^j5ip\p, s) = Aqs^, 
the fiavor-,S77(3) octet and singlet part, a 8 and a = AS are given by 

a 8 = Au + Ad- 2 As, AS = An + Ad + As. (3) 

The scale-dependent density AS evolves as 

Taking the difference between and T™ leads to the QCD Bjorken sum rule: 

If" = f dx\£(x, Q 2 ) - g?(x, Q 2 )] = \g A [l - + 0(a 2 s )]. (5) 

JO 7T 



The first order QCD correction was calculated in ITT], [L2| and the higher order correc 



tions were given in 13 



Now, the twist-4 operator gives rise to 0(1/Q 2 ) corrections to the RHS of (0). Their 
flavor decompositions are just the same as those of leading twists. From renormaliza- 
tion group analysis, their logQ 2 dependences take the following form in leading-log 
approximation. 



™ (n2 , 8 r 1 1 a s (Q 2 o) Y lNS, * PO . 1 (a.{Q 



(6) 



9Q 2 L L 12" 36 JOJ WW 9 JU \a s (Q 2 )J 

where /o, /3 and /§ are the twist-4 counter parts of ao, 03 and ag. 

(p, a |iT>, S )=/ i S<7 (7) 

/i's are scale dependent and here they are those at Qfo. Here we assume that there 
exists only one physical operator for each flavor as we show later. 7^ and 7^ are 
the coefficients of the one-loop anomalous dimensions for flavor non-singlet and singlet 
operators respectively. They are obtained from the renormalization constants of the 



corresponding operators. The magnitude of these corrections depends on the reduced 
matrix elements fi which are not calculable by perturbative QCD. 

On the other hand, target mass effects also give rise to power corrections to these 
sum rules. They can be estimated in full order of M 2 /Q 2 by taking the difference 
between the Nachtmann moments and the usual 1st moments fl4jl . 

ATi = Ti(M 2 ^ 0) - Ti(M 2 = 0) 

= L ^ { 9^ + 9 ^ V 1 + " Q } " 3 7o ^-— gM 2 m 

It should be noted that we need also g^x, Q 2 ) to estimate Ar 1; and this is because we 
could not obtain their values until the recent experiments of g^x, Q 2 ) came out. From 
these experiments [0, [| we get, 

Ar? = -0.001 ± 0.002, Ar™ = -O.OOl ± 0.002 (Q 2 = 2.5GeV 2 ). (9) 

So we can conclude that the target mass effects are negligible for the Bjorken sum rule 
and Ellis- Jaffe sum rule even in this lower Q 2 region because they amount to less than 
a few percent of IY 

In the following sections we calculate the anomalous dimension of the higher twist 
operators which determine the logarithmic Q 2 dependence of the twist-4 terms in the 
first moments of gi(x,Q 2 ). 

3 Flavor non- singlet part 

We now consider the renormalization of the operators. Let us consider the off-shell 
Green's function of twist-4 composite operators keeping the EOM operators as inde- 
pendent operators. Thus we can avoid the subtle infrared divergence which may appear 
in the on-shell amplitude with massless particle in the external lines. 

Another advantage to study the off-shell Green's function is that we can keep the 
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information on the operator mixing problem. And further, the calculation is much 
more straightforward than the one using the on-shell conditions. 

From general arguments it is known that there appear three types of operators 
which participate in renormalization of gauge invariant operators. B [15|. 



(1) Gauge invariant operators Ri which appear in the operator product expansion. 
We call them physical operators for they have non-zero physical matrix 
elements. 

(2) EOM operators E^ whose physical matrix elements vanish. 

(3) BRS-exact operators Ai whose physical matrix elements also vanish. 



These operators mix with each other through renormalization, and thief renormal- 
ization matrix is to be triangular. 

{ Ri,\ ( Zrr Z re Z R a \ f Ri \ 
Et = Z EE Z ea E t . (10) 

V A )r V o o z AA ) {a, ) b 

We should note that only Zrr has physical meaning among these renormalization 
matrix elements because of the triangularity. 

Now for renormalization of the twist-4 operators, at first we need to identify the 
physical operators and other operators which mix with them. As can be seen from 
the dimensional counting, there is no contribution from the four-fermi type twist-4 
operator to the first moment of gi(x, Q 2 ). The only relevant twist-4 operators turn out 
to be of the form bilinear in quark fields and linear in the gluon field strength. This is 
in contrast to the unpolarized case, where both types of twist-4 operators contribute. 

Operators which mix with each other need to have same properties such as dimen- 
sion, spin and other quantum numbers. The relevant operators in our case has the 
following properties; It is dimension 5 and spin 1. Its parity is odd and it has to satisfy 
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the charge conjugation invariance. The flavor non-singlet operators are bi-linear in 
fermion fields. We consider the gauge variant EOM operators as well, but BRS-exact 
operators don't mix because of flavor symmetry. 

The parity and charge conjugation condition requires the relevant operators to have 
an odd number of gamma matrices together with 75 or e^ upcr . The possible twist-4 
operators bilinear in tp and ip are of the form, O a = ipM a ip where dimOo- = 5, namely 
dimM, = 2, and we have M a = l5la D 2 ,gG^Y,l5(D a p+ pD a ), l5 (d a p+ pd a ). 
Hence we have the following operators which satisfy the above conditions: 

Ri = -^757 a £ 2 Vs 

R° 2 = g^G*^, 

E{ = ^75 pj a 0i> - ^ 75 £T p^ - ^ 75 pD a ij, (11) 

K = ^7s«9 CT p^J+^pd^, 

El = ? 757 a 9 W + ^75 P @l a ip, 

where = — igA®T a is the covariant derivative and G pv = \e ^vapG^ is the dual 
field strength. For example, charge conjugation (C) and parity (V) transformations 
read as follows: 

CR-C- 1 = -g4> T C-\-G^) T 1 ,£i r 
= -gViG^fllt 
= gW»lvti> 

= R2, (12) 
VRSV- 1 = #7o(-Gp 7 M7oV' 

= -#^7^ 

= -Rl- (13) 
Here one should note that not all of the above operators are independent, as in the 
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case of twist-3 operators ||16|| , and they are subject to the following constraint: 

R^R' + El (14) 

where we have used the identities, = ^{j^p} and [D^,D V ] = —igG^. Therefore 
any four operators out of flllD are independent and they may mix through renormal- 
ization. 

Now, Ei contains three gamma matrices, and when we rewrite the product in terms 
of one gamma matrices, there occurs a mixing from the gauge variant operator which is 
not an EOM operator. Here we take (R2, E\, E2, E$) to be the base of the independent 
operators. The only operator which really contribute to the physical matrix element 
responsible for 5V\ is R2. This twist-4 operator corresponds to the trace part of twist-3 
operator, (i^)^^ — Q^G^^^yij) — (traceterms), but there is no relation between 
the base for the twist-4 and that for the twist-3 operators. 

If we take this basis we have the following renormalization mixing matrix in the 
form of 



/ R2 






( Zu 


Z12 


Z\z 


Zu 


\ 


( R2 






Ei 







Z22 


Z23 


Z24 


Ei 




E 2 












Z33 







E 2 






V £3 


J 


R 


V 








Z44 


J 


V ^3 


J 


B 



(15) 



according to the general arguments in the following. (1) The counterterm for the EOM 
operators are given by the the EOM operators themselves. This is because the on-shell 
matrix elements vanish for the EOM operators ]T7| . (2) A certain type of operators do 
not get renormalized. And if we take those operators as one of the independent base, 
the calculation becomes much simpler. (3) The gauge variant operators also contribute 
to the mixing. 

Now we turn to the calculation of the renormalization matrix of this set of the 
operators. Since the R2 operator does not contribute to 2-point functions but to 3- 
point functions with quarks ?/>, ip and a gluon, A as external lines, at the tree level, we 
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only consider the following 3-point Green's function. 

rg A = (o\T^( P ')A;(im P )o a m\o) lpI , (i6) 

where these fields and the coupling constant involved are the bare quantities. Here we 
employ the dimensional regularization and take the minimal subtraction scheme. Note 
that we can not take all the external lines on-shell except for a singular configuration 
of the momenta. 

The Feynman rules are the following. 



/ tree 



= 9B1 5 1«(P + P) P T a - ig B e^l a l P T\ 

/ tree 

rlf) = -g B l,l P (p + P%T\ (17) 



tree 



Ilf ) = 9B%9a P ^+ f")T a - g B %1 P (p + P%T a 

/ tree 

-9Bl 5 la(p + p') P T a + ig B e apa(3 l a ^T a . 

The Feynman diagrams contributing to these Green's functions are in Fig.l. The 
one-loop radiative corrections lead to 

(rthi.ioop = {i + ^[-|c 2 (G)]}< toe wVr° 

+ -^-\ C ^ R ) + \c 2 (G)}g Bl5la (p + P') P T a 
e lb7T z 5 4 

e I6n 2 3 4 

1 9% 1 
e 167r 2 3 

It should be noted that the tensor structure %gB^a P a^ a l^T a is also included in (t^ a ^J 
and ( r^t A ) . So we can not connect directly the coefficient of the first line of (|T 

V 3 / tree 

with the renormalization constant of R 2 . We need to separate this coefficient into 
the i?2-part and the EOM-parts. In this sense EOM operators also contribute to the 

8 



+ -^[\c2{R)-- A C 2 {G)]g Bl <MP + P\T a (is. 



+ 777^o^(^75^+ P")T° 



tree 



physical quantity Jns- Rearranging the results of properly, we obtain; 

(rthi-ioop = {^ + l^[-lc2(R) + c 2 (G)]\(r^\ ree 



1 9b r , 1 



C 2 {R) + -C 2 {GW%^\ 



+ ^[-^ 2 (i?) + ^ 2 (G)](r|r)tree 
1 9% 1 ^ 



Now we introduce renormalization constants as follows: 

A b = Z l /2 A% = Zl /2 tp R , g B = Z g g R = Z x Z^ l2 g R . 

The composite operators are renormalized as 

(O i ) R = J2Z lJ (0,) B , 

3 

and the Green's functions of the composite operators with ip, ip and A as the 
lines are renormalized as follows: 

3 

For example, (T R2 ) R reads 

a/2 r„ r /. 3,5 



(r R2 ) fl = z 2 z 3 1/z ^ 11 (( v i + 7 (--c 2 ( J R) + c 2 (G))j(r /?2 ) tree 

+ |(-|c 2 (i2) + 1(^(0))^)^ 

+ -^ C 2(R) (r£ 3 ) tree J 

a a a 

+ jz 12 [T Entree + ( r £2Jtree + (r^sJtree 

where a = — Since we have in the Feynman gauge 
I67H 



Z 2 ZI' 2 = {l - ^[C 2 (i2) + C 2 (G)]} 
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the above equation (^3|) becomes 



a 



e 



(T R2 ) R = l + ^[ Zll -^C 2 (R)](T 9 ^ 



3 v /J V K2 /tree 



(Z12 - \c 2 {R) + -C 2 {G)) (rg) 



a 

+ f(^-^W + ic 2 (G))(r|») tree 

'tree 



+ + ic 2 (fl)) (rg) (25) 



which should be finite. z$j is defined as, 



■Z'tj = + ~ Zi ii (26) 



and (rf^) tree denotes (r fi2 ) tree with g B replaced by g^. For T^ A , there are additional 
diagrams due to tree quark- ant iquark vertices (Fig.2). 



(r^Jjj - ^2^3 1/2 [^22 (r^Js + ^23 {Fe 2 ) B + Z 24 E 3 , 



= i + ^h2-c 2 (i?)-ic 2 (G)](rg) tree 
+ f[,23 + 2C 2 (i?) + ic 2 (G)](rg) tree 

+ ^[z 24 -C 2 (R)- 3 -C 2 (G)](T%) tree . (27) 

Further, the EOM operators like E 2 and E3 which are of the form E = ipB-p=, where 

dip 

B is independent of fields, do not get renormalized |J. This can be seen as follows. 
Z 2 Z 3 1/2 (0\T U(x l )A;(x 2 )^(x 3 )^B 5 ^(y) ) j |0> 

= -iZ 2 Z 3 l/2 J V^VA^( Xl )A;(x 2 )^(x 3 )^B 6 -^(y) 

= -iZ 2 Z 3 1/2 J VipVipVA [ilj( Xl )A a p (x 2 )^(x 3 )Be iS 6\0) (28) 

+^x 1 )A a p (x 2 m y )e tS 6 4 (y-x 3 )} 
= -iZ 2 Z 3 1/2 (0\T U( Xl )A a (x 2 )^B(y)) \0)5 4 (y - x 3 ), 



where the last equality holds as <5 4 (0) vanishes in the sense of dimensional regulariza- 
tion. If B does not contain any field, the RHS of (|29|) becomes finite by the wave 
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5 S 

functional renormalization for ip,ip,A. Therefore E = ipB^ is a finite operator. 

dip 

Explicit calculations also indicate 

(r^ = z 2 z l 3 /2 z 33 {v E2 ) B = (1 + ^ 33 ) (rg) tree , (29) 

and 

(v Ea ) R = z 2 z] /2 z A , (v Es ) B = (1 + ^ 44 ) (rg) tree • (30) 

From the fmiteness of (0), (pTf), (p^), (|30"D, we get the following results for the 
renormalization constants: 

zii = |C 2 (i2), i zi 2 = |C 2 (i2) - iC 2 (G), 

^13 — 3 C2[R) — ^C 2 (G), = — gC 2 (i?), , V 

^22 = C 2 (i2) + |C 2 (G), z 23 = -2C2(i?) - ±C 2 (G), 

^24 = C 2 (i?) + fC 2 (G), Z 3 3 = ^44 = 0. 

This result is in agreement with the general theorem on the renormalization mixing 
matrix [§, || [15]. We should note that gauge variant EOM operator is also necessary 



to renormalize the physical operator. 

We now determine the anomalous dimension of R 2 operator. In physical matrix 
elements, the EOM operators do not contribute and we have 



(phys|(^) B |phys) = Zii 1 (phys|(i^) fl |phys) 

(phys|(^) R |phys). (32) 



9 2 1 



1 - -^--C 2 (R) 
16vr 2 e3 V J 

Therefore the anomalous dimension 7^5 turns out to be 



-f NS (g) = Zun—{Z l4 



d_ 



2z n + 0(g 4 ) 



9 . 0.9,, +D(r 4 

16vr 2 



16tt 
and 



,jl l NS + O(9% (33) 



X6 

= 2z n = -C 2 (R), (34) 
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which coincides with the result obtained by Shuryak and Veinshtein JTI| using the 
background field method. 

Including the twist-4 effect the Bjorken sum rule becomes 



dx 



in the case of QCD. 



^(*,Q 2 )-<0r,Q 2 ) 



71 



8 x /logQ 2 /logA 2 \~ 32/9/V 



r./; 



9Q 2 ^ \\ogQ 2 /\ogA\ 



(35) 



4 Flavor singlet part 

So far we have considered the flavor non-singlet part. Now we turn to the flavor 
singlet component. We should generally take account of gluon operators and BRS- 
exact operators as well in this case. 

At first we see whether there exists other physical operators in addition to R,2- The 
possible twist-4 and spin-1 operators are the following: 

G^frGii = G«D^G$ + G«D»G° = 0, 



G«Z>G^ = (from Bianchi identity), 



G aa D»G 



fia 



R3, 



(36) 



where D ab = d fl 5 ab + f abc T c . Further, i? 3 , a gluon EOM operator and a BRS-exact 
operator may enter into the mixing. 



E 4 = G aa D^G^ + -G^d a (dA)-gf abc GZ(d^ b u c )-g^G a ^, 

a 

A ee -G^d a (dA)-gf abc GT(d a ^c). (37) 
a 

Then we have a trivial relation among them, 



R 3 = Ea + R 2 — A. 



(3f 
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Therefore it happens that the physical operator is still only R 2 even in the flavor singlet 
case. This is characteristic to the lowest-spin case of higher-twist, in which only a few 
tensor structures are possible. 

It should be noted that we have only to consider the mixing between R% and 



E G = G a °D»G m -g^ a G a ^, 



(39) 



as long as we want to obtain the physical matrix element Zf x where we do not have to 



consider Green's functions with ghost operators in their external legs [18 



The mixing between R 2 and Eq can be studied by computing the Green's functions 
with two-gluon external lines shown in Fig.3. Then we have 



= ^<-! w '> x (^)w 

where Nf denotes the number of fermions. (r^e^ is given by 

( r ^) t ree = fix ~ QfiQx) ~ ie° XaP q*(q 2 9fi P ~ W P ) 

Now we introduce the renormalization constant Z15 as 



(40) 



(41) 



{R 2 ) R = Z n (R 2 ) B + Z 12 {E 1 ) B + Z 13 {E 2 ) B + Z 14 (E 3 ) B + Z 15 {E G ) B . (42) 



So we get 



i=l 



a 



N f) (r^) tree + z 15 (r^) tree 



e y 3 

From the finiteness of the above expression, we have 



On the other hand, the relation fl22|) becomes 

3 



Rp A ^j 



R 



7 7I/2 

Zj 2 Zj 3 



'11 



nr) B + E zu \r» A ) B + z 15 (vZ a ) b 



(43) 



(44) 



(45) 
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in this case. Since (Tf^ A )B is just the same as the flavor non-singlet case, we can easily 
extract Zf x from the result of Z^[ s . 



namely, 



3 

hence we obtain the exponent for the singlet part 



Z?i = Z^ S + 2 ~N f , (46) 



7° = 1°ns + t N f' ( 47 ) 



2(3 2(3 3 fa /3 V 9 3 
Again we reproduce the result of |TD|. Substituting these results into (0), we get 

™m C°»*-< /logOVlogA^\- <32/9Al+2 "' /3> ,^ 

sr a _ t = ( log< 3 S/logA2 J ■ m 

where C twlst ~ 4 denotes the coefficients including the reduced matrix elements, fa. 



5 Concluding remarks 



In this paper, we have investigated the renormalization of the twist-4 operators relevant 
for the Bjorken and Ellis- Jaffe sum rules. We obtained the anomalous dimensions 
by considering off-shell Green's functions in which the role of EOM operators was 
important. We have also seen that the renormalization matrix satisfies the general 
theorems for the renormalization of composite operators in the gauge theory. Since 
there appears only one physical operator, our case seems to be the simplest example 
which contains non-trivial structures. 

Phenomenologically whether the twist-4 effects has significant contributions de- 
pends on the value of their reduced matrix elements. Estimations of their values were 



given by various methods. For example, from the QCD sum rules [19|, STi amounts 



to about 4% of Ti n , and the chiral-bag model [^(J predicts similar value, but some 
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other estimations predict twice or much larger values pi| . Furthermore the reduced 
matrix elements of twist-4 operators are considered to be related to the renormalon 



ambiguity of the coefficient functions of the leading twist parts So we can not 

give any definite predictions about the reduced matrix elements until some consistent 
treatments are developed for these problems. On the other hand, the anomalous di- 
mensions of twist-4 operators are free from the renormalon ambiguity and have definite 
values. 

The numerical values of the exponents of log-terms in 5rf' n are so small that the 
log Q 2 scalings are quite hard to observe in experiments. Since the anomalous dimen- 
sions are considered to become large as the spin of the operators increases, the scaling 
of the twist-4 terms are probably easier to measure in some higher moments . However 
to obtain the anomalous dimensions in spin-n case by our method is not so easy be- 
cause many gauge variant EOM operators may enter into the operator mixing. In the 
twist-4 case, it is hard to simplify these structures by the method used in the twist-3 
23 , [Z4|] in which many tree amplitude of different EOM operators are projected 



case 



out to a single one by a null vector Q p . Since the twist-4 operators have a pair of 
Lorentz indices contracted, there still remain so many independent tree amplitudes af- 
ter the contraction by fl p and this complicates the situation. So we need some further 
developments in the techniques to disentangle the mixing problem. 
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